Robust utility maximization in a discontinuous filtration 



Monique Jeanblanc * Anis Matoussi ^ Armand Ngoupeyou 

University of Evry Val d'Essonne University of Maine University of Evry 

^ ■ January 20, 2012 

O 
(N 

■ Abstract 



We study a problem of utility maxunization under model uncertainty with infor- 
mation including jumps. We prove first that the value process of the robust stochas- 
tic control problem is described by the solution of a quadratic-exponential backward 
stochastic differential equation with jumps. Then, we establish a dynamic maximum 
principle for the optimal control of the maximization problem. The characterization 
of the optimal model and the optimal control (consumption-investment) is given via a 
• forward-backward system which generalizes the result of Duffie and Skiadas [14] and 

I El Karoui et al. [18] in the case of maximization of recursive utilities including model 

with jumps. 
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CN ! 1 Introduction 

The utility maximization is a basic problem in mathematical finance. There is a 
^ I huge literature on this problem based on two approaches: the Bellman approach and 

■ the martingale one. Karatzas et al. [27] study a consumption-investment problem in a 

more general case. Their model takes into consideration the inherent non-negativity of 
consumption and considers bankruptcy constraint. When the risky assets are modeled by 
geometric Brownian motions, the value function is determined explicitly as well as the 
optimal consumption and the investment strategy by solving the Bellman Equation and 
using a verification theorem. Martingale methods were introduced by Karatzas et al. [28] 
and Cox and Huang [10]. They characterized the optimal consumption- portfolio policies 
when there are non-negativity constraints on both consumption and final wealth. Unlike 
the nonlinear Bellman equation, they give a verification theorem which involves a linear 
partial differential equation. In all these works, the authors assume that the underlying 
model is exactly known. 
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Our paper deals with the problem of utility maximization from a terminal value and 
an intermediate control under model uncertainty. Uncertainty refers to the case in which 
a decision marker does not know the probability distribution governing the stochastic 
nature of the problem she/he is facing. This uncertainty is captured by using capacities or 
sets of probability measures over the space of state of the world. The set of such probability 
measures on some measurable space is called by economics, objectively rational 

beliefs, and each element of such set is called beliefs on T that the decision maker is able 
to justify on the basis of the available information. The incompletness of information in 
then captured by the fact that we consider a non singleton set of probability measures. 

In the mathematical finance literature there are two approaches to solve robust util- 
ity maximization problems. The first one relies on duality methods such as Quenez [39], 
Gundel [20], Shied and Wu [41], and Shied [40]. The second approach, which is followed 
in this paper, is based on the penalization method and the minimization is taken over all 
possible models as in Anderson, Hansen and Sargent [2] and Hansan et al. [21]. These 
authors have introduced and discussed the basic problem of robust utility maximization 
when model uncertainty Q is penalized by a relative entropy term with respect to a given 
reference probability measure P. Both papers are cast in a Markovian setting and use 
mainly formal manipulations of Hamilton-Jacobi-Bellman (HJB) equations to provide in- 
sights about the optimal investment behaviour in these situations. Moreover, Skiadas [44] 
follows the same point of view and gives the dynamics of the control problem via BSDE 
in the Markovian context. More precisely, Skiadas [44] points out that the BSDE coincides 
with the one describing a stochastic differential utility; hence working with a standard ex- 
pected utility under (a particular form of) model uncertainty is equivalent to working with 
a corresponding stochastic differential utility under a fixed model, see also Duffie and Ep- 
stein [13] and Duffie and Skiadas [14] for more about stochastic differential utilities. We 
have to mention the interesting works of Maccheroni, Marinacci and Rustichini [33, 34] 
which have studied preferences and robustness by using variational technics (principle). 

More recently, Bordigoni, Matoussi and Schweizer [7] have studied this robust prob- 
lem in more general setting (non Markovian approach) by using stochastic control tech- 
nics. They considered the robust maximization problem as following : 

sup inf U((V',c),Q) (1) 

7r,c QGQ 

where runs through a set of random variables, c through a set of control processes and 
Q through a set of models (measures), and where the criteria U ((V', c) , Q) is the sum of 
a Q-expected utility and (a penalization term) a term relative to the relative entropy. They 
have solved only the minimization problem and they proved the existence of a unique 
Q* optimal minimizer model. Moreover, in the case of continuous filtration, they use the 
dynamic programming Bellman principle to show that the value function of the stochastic 
control problem is the unique solution of some generalized BSDE with a quadratic driver. 

Bordigoni [6] studied partially the maximization problem by using classical opti- 
mization arguments and assuming some conjectures. She derived the Gateaux differential 
of the Lagrangian associated to the optimization problem. She has obtained necessary and 
sufficient conditions that must be fulfilled for the optimal strategy in a complete market 
in the case of consumption/investment problem. 

More recently, Faidi, Matoussi and Mnif [15] have studied the maximization part of 
problem (1) by using the BSDE approach as in Duffie and Skiadas [14] and El Karoui, Peng 
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and Quenez [18] in the case of continuous filtration. We mention also that there is another 
approach based on Ekeland Variational principle to obtain a dynamic maximum principle 
for recursive utility optimization problem (see Ji and Zhou [24, 25]). 

Quadratic BSDE's in the continuous framework have been introduced by Kobylan- 
ski [29] and then generalized in [31, 8, 9, 35]. The extension to quadratic BSDEs with jumps 
were introduced by Becherer [5] and Morlais [36, 37] in the case of bounded terminal con- 
dition. In our paper, we prove first that in the case of discontinuous filtration (information 
including jumps), the value process V of our stochastic minimization problem in (1) is 
described by a particular class of quadratic generalized BSDE with jumps. Then, in or- 
der to tackle the maximization problem, we prove a comparison theorem for this class 
of quadratic BSDE with unbounded final condition, which allows us to derive a maxi- 
mum principle. More recently. El Karoui, Matoussi and Ngoupeyou [19] have obtained an 
existence result of a general class of quadratic BSDE's with jumps and unbounded final 
condition. 

Moreover, we give a necessary and sufficient condition of optimality and a charac- 
terization of the optimal control {if^* , c* , Q*) as the solution of a forward-backward system. 
Finally, we prove in the case of logarithmic utility function of the control process c that 
the optimal control (0, c* ,Q*) is characterized via the unique solution of some forward- 
backward stochastic differential equation. 

By using BSDE technics, we generalize the characterization of optimality obtained by El 
Karoui, Quenez and Peng [18] in the framework of robust case and including model with 
jumps. Our results may also be considered as a generalization of the works of ([14, 42, 43]). 

The paper is organized as follows. In Section 2, we describe the model and the 
stochastic control problem. In Section 3, we characterize the optimal model measure for 
the minimization problem via a solution of quadratic BSDEs. We prove a comparison and 
a concavity results of the solution of our BSDE with respect to control parameters. In 
Section 4, we derive the necessary and sufficient conditions that must satisfy the optimal 
control and then we establish the dynamic maximum principle which characterizes im- 
plicitly the optimal strategy (V'* , c* , Q* ) as a solution of a forward-backward system. For 
a specific choice of utility functions, the value function is given in Section 5 in terms of 
the optimal plan as a solution of a forward-backward system. Finally, Section 6 contains 
a technical proof concerning a regularity result of our quadratic-exponential backward 
stochastic differential equation. 

2 The IVIodel and the stochastic control problem 

We consider a filtered probability space (0, Q, G, P). All the processes are G-adapted, 
and defined on the time interval [0, T] where T is the finite horizon. We recall that any 
special G-semimartingale Y admits a canonical decomposition Y = Yq + A + M^''^ + M^'*^ 
where ^ is a predictable finite variation process, M^''^ is a continuous martingale and M^''^ 
is a discontinuous martingale. 

Assumption A 1. We make the following assumptions: 

1. For each i = 1, . . . ,d, W is a counting process and there exists a positive adapted process 
\\ called the P intensity of H^, such that the process with 




(2) 
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is a martingale. We assume that the processes H\i = 1, . . . ,d have no common jumps. 

2. Any discontinuous martingale M'^ admits a representation of the form dMf = X^f^^ yldNJ: 
where y\i = 1, . . . ,d are predictable processes. 

This hypothesis is satisfied in the case where the filtration is generated by a continuous 
martingale and an independent inhomogeneous Poisson process, or under regularity as- 
sumptions in the case where the counting processes are stopped after the first jump, as it 
is done in credit risk, as in the following example: 

Example 1. (see [23]) For each i = 1, ... ^d, let HI = l{T-.<t}, where Ti is a random time and 
F(rj = Tj) = Oj-j 7^ j. Let G be the smallest right-continuous filtration which contains the 
filtration ¥^ generated by a p dimensional Brownian motion B and the filtration generated by 
the processes H\ Under the assumption that P(rj G dOi; i = 1, . . . , d\Fl^) ~ P(ri G dOi; i = 
1, . . . , d), then any local G-martingale ( = {Ct)t>o admits the following decomposition: F-a.s, 

Ct = Co + /* Zs ■ dWs + fus-dN, > (3) 

JO JQ 

where W is the martingale part of the G-semimartingale B , Z := {Z^, . . . , Z^) and U := 
(U^, . . . , V^) are G predictable processes. Furthermore, if C, is square integrable 

< OO. 



E[ / iZ.pds] < OO, Ve[ /" \Ui\^Xids\ 
Jo Jo 



We denote by |X| the Euclidean norm of a vector or a row vector X. We give now some 
notations and definitions : 

Definition 1. 

L^^P is the space of all Qx-measurable random variables X with ^ [exp (7jX|)] < oo, V7 > 0. 
is the space of progressively measurable processes X = {Xt)^^rp with 

^ [exp (7 esssupQ<j<2-|Xj|)] < 00, V7 > 0. 
is the space of progressively measurable processes X = {Xt)^^rp such that 



E^ 



T 



exp I 7 / |ds 







< 00, V7 > 0. 



7V1^(P) is the space of¥-martingales M = {Mt)^^r^ with Mq = and [ supo<t<T \Mt \p] < 00. 

d I.J' 



£^(A, P) is the space ofW^-valued predictable processes X such that E"^ / {Xl)'^\\ds 
V^iF) is the space ofR-valued predictable processes X such that WF" ^ X'^ds 

5^(P) is the space of all R-valued predictable processes X such that E' 
Mq locO^) is the set of continuous local martingales null at zero. 



< 00. 



sup \X. 

0<t<T 



< 00 

2 



< 00. 
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Definition 2. 

For any probability measure Q on ($7, Gt), 

[ +00 otherwise 

is the relative entropy o/Q with respect to P. We denote by Qf (resp. Q^j the space of all probability 
measures Q on (f], Qt) with Q P (resp. equivalent to Pj on Qt and i?(Q|P) < +oo. Note that 
the reference probability measure P belongs to Qj. 

2.1 The robust optimization problem 

We define a discounting process 5"^^ := e~ ^"'^^ for all t € [0, T] where 5 is a non- 
negative adapted process. For Q € Q/, we denote by Z'^ = (^?)o<t<T cadlag P- 
martingale) its Radon-Nikodym density with respect to P. For Q G Q/, we denote by 
o<t<T Radon-Nikodym density with respect to P. Let C/ be a given process 
(the cost process) and IJt a given random variable (the terminal target). The robust utility 
maximization problem 'P{U, Ut, /3) is to find the infimum of r(Q) over the set Qf where 



/ S^Usds + Sl^UT +/3E^ / 6sSilnZ^ds + Sl^lnZ^ 
Jo J Uo J (4) 



■0,T 



The first term in the right-hand side of (4) will be linked, in the following section, to the Q- 
expected discounted utility from target and cost process. The second term is a discounted 
relative entropy term and /3 > is a given positive constant which determines the strength 
of this penalty term. Note that the optimal probability Q for the problem V{U, Ut,P) 
is optimal for the minimization problem 7^(11^ ,u!^,l) where = U/f3,tJj, = Ut/P, 
therefore, we shall restrict our attention to the problem V{U^ Ut) ■= V{U, Ut, 1). 

Assumption A 2. For a more precise formulation of our problem, we make the following further 
assumptions: 

i) the discount rate 6 is a non-negative bounded process, more precisely there exist e > such that 
for any t > 0, < e < (5j < ||5||oo, a.s. where \ \d\\oo is a constant. 

ii) the cost process U belongs to D°^^ and the terminal target Ut is in L'^^p. 
Hi) the process := A^ds is uniformly bounded, i.e., < C, a.s.. 

Remark 1. The assumption U belongs to Dj^^ implies that E'^ [ \ Us\ds] < oo for allQ e Qf. 
Indeed, applying the following estimate: 

xy < ylny — y + e'^ , for allx gM., y > 0, 

we get 

f \Us\ds\ = E^[Z^ r \Us\ds\ <H{Q\¥)-E^[Z^] +E^[exp( f \Us\ds) 
Jo Jo Jo 
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Remark 2. The assumption Hi) is a technical hypothesis needed only in the proof of Theorem 4. 

We recall the existence result of the optimal probability measure for the minimiza- 
tion problem infQ^g^ r(Q) which was given in Theorem 9 and Theorem 12 in Bordigoni, 
Matoussi et Schweizer [7]: 

Proposition 1. Under Assumptions A1-A2, there exists a unique Q* ivhich minimizes T{Q) over 
all QeQf: 

r(Q*) = inf r(Q) (5) 
Furthermore, Q* is equivalent to F, i.e Q* G Q'j. 



3 The Optimal Model Measure and BSDE 

3.1 A BSDE description of the value process 

We use stochastic control techniques to describe the dynamic of the value process 
V associated to our robust optimization problem, via BSDEs. In a markovian framework 
(see Skiadas [44]) or in a semimartingale setting (see Bordigoni, Matoussi and Schweizer 
[7]), the authors have established that V is the unique solution of a backward stochastic 
differential equation (BSDE) with a quadratic driver. In our paper, the BSDE associated to 
our control problem is defined as follows : 

Definition 3. A triple of processes {Y,M'^''^,y) such that Y is a F-semimartingale, M^''^ is a 
locally square-integrable continuous local F-martingale null at and y = {y^ , ■ ■ ■ ,y^) an M*^- 
valued predictable locally hounded process, is called solution of quadratic-exponential BSDEs, if, it 
satisfies: 



(6) 



a ^ a 

dYt=[Y, g{yl)K -Ut + StYt] dt + -(i(M^'<=), + dM^^'' + ^ yldNl 

i=l i=l 

Yt = Ut 

where g is the convex function g{x) = + x — 1. Note that Y is a special F-semimartingale. 

Remark 3. In the case the filtration G is generated by some given brozvnian motion W and the 
jump process N , then the BSDE takes the form of standard quadratic BSDE : 

d ^ d 

dYt=[Y, aivDK -Ut + 5tYi^ dt + -\Zt\^dt + Zt.dWt + yldNi 

1=1 i=l 

Yt = Ut 

Such BSDEs have been studied recently and typically appear in problems from pricing-hedging 
derivative options by indifference pricing or/and maximization of expected exponential utility in- 
cluding jumps on the wealth portfolio; see for instance Becherer [5], Morlais [36, 37] , Lim and 
Quenez [32], Ankirchner, Blanchet-Scalliet and Eyraud-Loisel [1] and Schroder and Skiadas [43[ 
for some recent references. However, all existence and comparison results for such equations assume 
that the terminal value Yt is bounded; here, we relax this condition. In a recent work. El Karoui, 
Matoussi and Ngoupeyou [19] have obtained an existence result of a general class of quadratic 
BSDEs with jumps and unbounded final condition. 
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We first establish a recursion relation for solutions of (6) which implies the unique- 
ness of the solution : 

Proposition 2. Let {Y,M^''',y) G Dl""^ x M^qi^^(F) x C^iX,]?) be a solution of the BSDE (6). 
Then, Y satisfies the following recursion equality: for any stopping time r valued in [t, T], 



Y 



•In 



exp 



Yr + / {6sYs - Us)ds 



(7) 



Moreover the BSDE (6) admits at most one solution which belongs to x Mq 



loc\ 



Proof: (i) Assume that {Y, M^^", y) is a solution of (6), and define Xt = Yt-YQ- fl{6sYs 



Us)ds and Zt = e"^'. Ito's formula leads to dZt = Zt- \ -dM^''' + Y!i=i (e^*'* - ij dNl 
Hence, Z is a non-negative local martingale. Assuming that Z is a martingale, one obtains, 
for t < r < T: 



-Yt 



E^ 



exp {-Yr + ^ {5sYs - Us)ds^ 



(8) 



In general, we use a localizing sequence Tn in order to have the P-martingale property and 
thus obtain (8) with Tn At instead of r. Then by the integrability Assumption 2 and the 
assumption that Y E D^^^, we obtain a P-integrable upper bound for the right-hand side 
of (8) and letting ?i go to infinity, by dominated convergence we obtain (7) for r. 
(ii) Uniqueness of the BSDE (6); Assume that {Y, M^'^, y) and {Y, M^^'^, y) are two solutions 
of (6) in D^f^^ X Ml^i^^{¥) x C^{\, P). Suppose that, for some t G [0, T], the set A = {Yt > 
Yt} G Qt satisfies P(A) > and define r = inf{s > t\Ys > Ys}, so that Yr > Yr. Since 
Yx = Yp, one has t <T, and: 



{5sYs-Us)ds-Yr> / {5sYs-Us)ds-Yr on A, 



then from the recursion relation (7), it follows that 



exp {-Yt) = W 



exp 



5sY, - Us)ds - Yr 



> exp {—Yt) on A 



which implies that Yt <Yt on Ain contradiction with the definition of A; therefore Y and 
Y are indistinguishable. □ 

Remark 4. In the case (5 = 0, the process Y, part of the solution of (6) is given in a closed form as 



Yt 



InE^ 



exp —Ut 



U^ds 



In the case U = 0, we recognize the dynamic entropic risk measure studied, in particidar, by 
Barrieu and EL Karoui [4]. 

The main result of this section gives the BSDE description of the value process of 
our robust control problem in the case of information including jumps. It extends earlier 
works by Skiadas [44] (Theorem 5, pp. 482) and Bordigoni, Matoussi and Schweizer [7] 
(Theorem 12) (see also Lazrak and Quenez [30] and Schroder and Skiadas [42]). 
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Theorem 1. Assume (Al) and (A2). Then there exists a unique triple {Y,M^''^,y) e D^^^ x 
■Mq(P) X £^(A,P) solution of (6). Furthermore, the optimal measure Q* solution of (5) admits 
the Radon-Nikodyjn density Z^* = £{L) w.r.t. P where 

d 

dLt = -dUl'" + ^[e-y'^ -l^dNl Lq = 0. (9) 

i=l 

Proof: We divide the proof in three steps. We first prove that the value process V of our 
control problem is a P-special semimartingale, i.e. V = Vq + + with = [M^''^ + 
XliLi vldNf. Secondly, we prove that {V, M^''^, v) is solution of the BSDE (6). Finally, we 

show that {V, M^'^ v) € i^o^^ x yW|;(P) x ^^{X, P). 

Step 1: We embed the minimization of r(Q) in a stochastic control problem and we use 
mainly the martingale optimality principle from El Karoui [16] (Theorem 1.15, Theorem 
1.17 and Theorem 1.21) to get our result. To that end, we introduce a few more notations. 
Let I denote the set of all G-stopping times r with values in [0, T] . We define the minimal 
conditional cost 

J(t, Q) := Q - ess miQ,^T>iQ,r)^{r, Q') 



with r(T,0') := Er 



^OT + ^OT(Q')l^r andP(Q,r) = {Z^' \Q' G Qf and 



on Qr}- So, we can write our minimization problem as 

inf r(Q) =E^[J(0,Q)] 



■f 



by using the fact that Q = P on for every Q G Q/. A measure Q € Q/ is called optimal 
if it minimizes Q iH> r(Q) over Q € Qf. We know from Proposition 1 that there exists 
an optimal Q* which belongs to Qp hence, w.l.o.g., we restrict our attention to minimize 
over Q € Q J . For each Q € Q j and t ^ I, we define with obvious notation 



y(r, Q) := Q - ess mfQ,^j,^Q^,)EQ, l^U^^^ + ^r,r(Q') I Sr 

which is the value of the control problem started at time r instead of and assuming 
one has used the model Q up to time r. By using the Bayes formula and the definition of 
t(Q')' one can easily prove that V{t, Q) = V{t) does not depend on Q e Q'j. Moreover, 
comparing the definitions of V{t) and J(r, Q) yields for Q € Qj 

J(r, Q) = S^^V{t) + r S^Usds + f 5,5f In Zfds + In 
Jo Jo 

because we can also take the ess inf for J(r, Q) under P Q. From the martingale 
optimality principle proved in [7] (Proposition 13 pp.140), there exists an adapted RCLL 
process J"^ = {j'j^)o<t<T which is a right closed Q-submartingale such that = J{t, Q). 
Thus we can choose an adapted RCLL process V = {Vt)o<t<T such that Vr = V{t) = 
Vir, Q), P - a.s. for t e I and Q € Q}, and then we get for each Q G Q^, 

= 5V + / SiUsds + / 6sSi In Zfds + In Z^. (10) 
Jo Jo 

As P G and is a P-submartingale (from Proposition 13 pp. 140 in [7]), equation (10) 
yields that = S^V + J SgUgds. Thus F is a P-special semimartingale, i.e. its canonical 
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decomposition can be written as 

V = Vo + M^ + A^. 

Since is uniformly bounded from below and is a P-submartingale, hence is a 
true P-martingale and from (Al) = M^''^ + Yli=i v^dN^ where M^'*^ is a continuous 
F-martingale. 

Step 2: We now prove that {V, M^''^, v) is solution of the BSDE (6) where v := {v^, - ■ ■ , v'^). 

For Q € Q^j, we denote by L'^ the stochastic logarithm of Z'^, i.e., the P-local martin- 
gale such that dZ^ = Z^^dL^. From Assumption 1, the local martingale admits the 
decomposition dL^ = dL^''^ + Ylt=i ^\dNl, where L'^''^ is a continuous P-local martingale, 
and are predictable processes, obviously greater than — 1 and one has 

^ d d 

din Zf = dLf'' - -diL^^t + E 1^(1 + ^"t^dNi + 5^(In(l + £1) - l^Xldt. (11) 

i=l 1=1 

Using integration by parts formula, we obtain after some simple computations and using 
(10) and (11): 

dJ^ = Sf (^i-6tVt + Ut)dt + {dVt + dlnZ^)^ 

= Sf [i-5tVt + Ut)dt + dM^'" + dA'Y + dlf'^ - ^d{L'^'')^ 

d d 

+ EK + in(i + e^))dNi + J](in(i + e,) - i^xidt 

i=l 1=1 

From Girsanov's theorem, the processes (ND^^^ and (Aff )j>q defined as dNl = dNl — 
iiXidt, and dM^r = d^M^'" + hf") - diM^'" + L^'^ L'^''')t are Q-local martingales, hence: 

dJ^ = sf[i-6tVt + Ut)dt + dM^ + dAY + d{M^''^ + L'i''^,L'i''^)t-^d{L'i''), 

d d 

-I- EK + + mdNl + J2 i^M - 1) + (1 + ii) ln(l + Q) Kdt . 

i=l i=l 

In order that the process J^^ is a Q-submartingale for each Q € Q'j, we impose that its 
finite variation part is a non-decreasing process. 

AY = -ess inf [\us - SsVs)ds + (M^-'^ + L^'^ L^'^)^ - ]-{L^'')t 
2/ JO ^ 

+ E [\iiivl-i) + {i + ei)H^ + £l))Kds. (12) 
i=i 

To find the ess inf, we divide (12) in two parts, the continuous part and the discontinuous 
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part; hence we have two optimization problems: 

AY = [\6sVs-Us)ds-essmf{{M'''^,L'^'^)t + l{L'^^')t} 
Jo 2/ ^ 

- essi;?/E r(4K-l) + (l+4)ln(l + 4))Aids. 
2/ Jo 



It is proved in [7] that the first infimum is obtained for L^'"" = -M^'"" and -ess inf { (M^'^, L^''^) + 



1 ^2^Q,c^| _ l^jv/^'*^). The second part of the optimisation problem reduces to find the op- 
timal t, solution of ess inf (f * — 1) + (1 + il) ln(l + )) which is an easy task, the solu- 
tion being t/ = e-< - I, which leads to -ess inf - 1) + (1 + 4) ln(l + 
= e~^'> + vl — 1 = g{v\) where g{x) = + x — 1. Therefore, 

AY = f{5sVs-Us)ds + \{M''nt+ fY.g{v\)Kds. 
Jo ^ Jo 

It follows that {V, M^''^, v) is a solution of 

dVt = (stVt -Ut + Y, s{vi)Xi] dt + \d{M^''')t + dM^'" + vldNi Vt = Ut (13) 

\ i=l J i=l 

Hence, the optimal probability measure Q* is characterized by its Radon-Nikodym density 

d 

dZf' = Z^:dLt, dLt = -dM^'' + - l) dNl. 

1=1 

Therefore there exists a solution of the quadratic BSDE (6) . 

Step 3: In this step we prove that the solution {Y,M'^''^,y) of the BSDE (6) belongs to 
the required spaces. 

From Lemma 19 and Proposition 20 in [7], we have that Y belongs to D^^^. Let now study 
the space of the process M^''^. We consider the F-martingale: 

rT 



Kt := 



exp(^^ {5sY,-Us)ds-tJT^ 



Qt 



Using the fact that Y G D^^^, we obtain that the process K belongs to JW{f). Now, the 
recursive property leads to Kt = exp {^—Yt + /q((5sYs — Us)ds^ and it is not difficult to 
show that, from Ito's formula and the canonical decomposition of Y , 

dMj^'' = . (14) 

Kt- 

From Assumption 1, there exists A;* and M^''^ such that the martingale K equals Kt = 
Ko + Mf + Y.i=i /o ^idNl. Hence, from (14): 

(M^'^)t< r ^d{K^)t < {K^)t sup ^ 

Jo K^ 0<t<T K^ 

<(K^)rexp(2 sup jy^Kl + ||5||ooT) + 2 \Us\ds] 

\ 0<t<T Jo J 
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By BDG's inequalities, there exists a constant C such that for every p € [1, +00): 

rT 



{K'h + [ 
Jo 



< CE^i sup \Ktf) 

0<t<T 



(15) 



Since K £ M^iF), we conclude that M^'^ lies in the space Mq{F) for every p G [1, +00). 
We conclude, using again BDG's inequalities. Finally, let now characterize the space of 
jump process y. Using the recursive relation and the decomposition of the process K we 

get ln{Kt- + ki) - hi{Kt- ) = -yl hence: 



{e-yt - 1) dHi 



E' 



Since sup f-^) € L^(P) for any p € [l,+oo], using (15) and Cauchy inequalities, we 



0<t<T ^J^t 
conclude: 



In particular 



(e-^* - 1) dHl 



< 00. 



E^ 



{e~y^ - 1) Xldt 



E^ 



J4_ 



dHl 



<E^ 



sup T72 

.0<t<T / JQ 



{k\)dH, 



(16) 



< 00 . 



By using similar arguments, one proves that: 



E^ 



T 



(e^' - 1) Ajdt 



< 00. 



(17) 



(18) 



Moreover, by using the inequality |yp < 2(|e~^ - Ip + |e^ - Ip), My € M and (17)-(18) we 
conclude that the process y is in £^(A, P). □ 

Remark 5. The martingale part of the BSDE solution, i.e., M = -M^^" + J^f^^ /o(e"^* - l)dNi 
belongs to Ml{¥)for any p e [1, +00). Indeed, since M^'" G A^o(P), and (16) 



E' 



d 



i=l -^0 



< 00, 



and using BDG inequality, we obtain E^ (supo<j<T |Mj |p) < 00. 



3.2 Properties of the value process 

In this part, we establish a comparison theorem for our class of quadratic BSDE (6) which 
is a key point to derive the dynamic maximum principle for the maximization problem. 
We shall study the properties of the value process for a given pair {U,Ut)- 

Definition 4. For two random variables X and Y, we write X < Y for X < Y a.s. For two 

processes A and B, we write A < B for At < Bt^Mt G [0, T], a.s.. We write {X, A) < {Y, B) if 
X <Y and A< B. 
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Theorem 2. Assume that for k = 1, 2, {Y^, M''''^, y'^) is the solution of the BSDE (6) associated 
with (U^, U^). We denote — yi _ y2^ jji2 ._ jji _ jji fjvi ._ jji _ jj2_ jj^^^^ 



(19) 



where Q*'^ is the solution ofV{U'^, U^), i.e., the probability measure equivalent to P with Radon 
Nikodym density Z'^*' given by 



\ 1=1 / 



t- 



(20) 



In particular, if{U^,U^) < {U^, U^), one obtains < Y^, dP ® dt-a.e. 
Proof: We denote y*'^^ — yi,i _ yi,2 ^ji2,c ^ ^i,c _ ^2,c^ ^]^^^ ^^^^ 

Jt .^1 Jt .^^ Jt 



+ 



T 



T 



dM. 



12,c 



Note that, since M^''^ are continuous martingales. 



(21) 



Using the fact that the process {M^'^''^) is increasing and that the function g is convex we 
get: 

Jt j^]^ Jt Jt 



12,c 



dM' 



i=l 



Let N* and M*''^ be the Q*'^-martingales obtained by Girsanov's transformation from N 
and Af ^2,c^ where dQ*'^ = Z^*'^d¥ and where Z^*'" is given by (20), then: 



d 



1=1 -^^ 



y^^^dNi*- dUr 



rl2 



which impUes that < ^^*'^\ f^ e' ^^'^^'Ul'^ds + e' ^-d'Vr 
(^7^ U^) < i7|), then Y^^ < Y^ dF ® dt-a.e. 

We now prove a standard a priori estimates for the solution of BSDE (6) 



. In particular, if 

□ 
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Proposition 3. (A priori estimates) Let M'^'^, y'') be the solution associated with {U^, U^) 
for k = 1,2 where we assume that {U^, tJ^) < (f/^, U^), then there exists a constant C > such 
that: 



sup + + E 

o<t<r 



1=1 



IVt I K dt 



12|2 



(22) 

where A*'* = A*e~^''^ is the intensity process of W under the probability Q*'^. In the case 
{U'^, U^) < {U^,iJ^), the same inequality holds with Q*'^. 

Proof: Using Ito's formula: 



12\2 



12 



+ 2Y, 



12 



d 



+ d{W''^)t 



where dmart/ = 2Y, 



12 



1=1 
12,c 



dt + ^{yi^'^^fXidt + dmartt 



i=l 



dM^''"' + Yli=i vi'^^dNi^ +T,Li {yt^^y dNi corresponds to a local 
martingale. Assuming {U^, tj^) < (f/^, U^)' follows from the comparison Theorem 2 
that < y^. Using the relation (21) and the convexity property of the function g, we get 



Jt Jt 

i=i Jt 
Jt ■ 1 Jt Jt 



i=l 

Hence we finally obtain the following inequality: 

iY,''f + fd{M''ns < {UPf + 2 fYP [-6sY}^ + Ul^] ds -j^ f {v^^^Kds 
Jt Jt j^^^ Jt 

+ I dmartl 



t 

i 2 

where mart* is a Q*'^ local martingale and A" := A* e~^«' is the intensity of W under Q*'^. 
From the obvious inequality 

(Yt'r - j^{U}^)Y,'' > ^{U^'f 
and the non-negativity of 5, we deduce easily that 



-Y,''{S,Y,^'-Ur^)<^^iUPf (23) 
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Plotting relation (23) and using the fact that the process 5 is bounded below, there exists a 
constant C > such that: 



sup |y,12|2^^^^12,c^^^'^ f 
tG[0,T] -^-^ Jo 









f'^ I 12,i|2,i,*, 




Jo 







Permuting Y and Y and assuming {U , U^) > {U , Uf.) leads to the kind of inequality. 

□ 

As a direct consequence of the comparison theorem, we prove the concavity property for 
the BSDE solution. 

Theorem 3. (Concavity property) Define the map F : D'^^ x L°^p — > D'^^ as 

F{U, U) = V 

where {V,M^''^,v) is the solution associated with {U,JJ). Then F is concave, namely, for all 
9 G (0, 1) and (C/\ U^), (U^, U^) G ^ x L^^p ; 

F [eu^ + (1 - eu^ + (i - e)u^) > eF{u^, u^) + (i - e)F{u^, u^). 

Proof: Let (V'', M'''", v^) be the solution of BSDE (6) associated with ([/^ C7|) e L>f ^ x 
L^^P, then for any 6* € (0, 1): 

d{ev,^ - (1 - e)v^) = [6t{evt' + (i - e)v,') - {eu} + (i - e)uh] dt 

+ 9d{M^'^)^ + (1 - 9)d{M^^'')^ + d{9Ml''' + (1 - 9)Mt''') 
d d 

+ [evY + (1 - e)vY] dNi + ^ [e^K''^) + (1 - e)g{vY)\Kdt 



1=1 



1=1 



We recall the following general result: Let X and Y be two continuous martingales. Then, 
for all e G (0, 1), 9{X) + (1 - e){Y) - {9X + (1 - e)Y) is an increasing process. Indeed, we 
have: 

{9X + (1 - 0)Y) - 0{X) - (1 - 9){Y) 
= {9^ - 9) {X) + [{l-9f -{I- 9)) (y) + 29{l - e){X, Y) 

= e{9-i)[{x) + {Y)-2{x,Y)] =e{e-i){x-Y) 

Therefore using the convexity property of the function g we get: 

9Vl + {l-e)V^ < [eu}, + {l-9)U^)- j\5,{9V} + {l-e)V^)-{9Ul + {l-9)Ul)\ds 

- [ d{9M^'^ + (1 - 0)m2'^)^ - i diOM}^'' + (1 - 9)M}'') 
Jt Jt 

- E / (^^^'' + (1 - ^ys')dNl -Y^l 9{0vY + (1 - e)vY)\\ds (24) 
1=1 -^^ i=i -^^ 
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Let (F^ M^'^ v'^) be the solution of the BSDE associated with {OU^ + (1 - e)U^, OU^ +(1 - 
e)U'^) and set M^'^'^ = OM^^'' + (1 - 6')M2'^ and for i = 1, • • • , d, 9^'* = Ov^^' + (1 - 6')?^^'% 
then using (24): 

9v,^ + (i-0)v,2_^^e< r Ss{v! -{ev} + {i-e)v^))ds- r d{M''^^^<^)^+ r d{M'^^)^ 

Jt Jt Jt 



Jt .^^ .^^ 



9f'^-z;f'^)diV: 



Using (21) and the convexity property of the function g we get: 

T 



ds 



+ E 



i=l 



0,i 



[ d(Mr'^'^ - M^) - / {vt^^ 

Jt .^1 Jt 



^}^'*-t;r)diVl 



Therefore, we have: 



ev,' + (1 - e)v,^ - < j"^ [5s{v! - {ev} + (i - e)vj)\ ds 



Let Q*'^ be the probability measure equivalent to P with Radon-Nikodym density dZ^ = 

—dM^''^ + Yli=ii^~~^''' ~ ^)dNl , then using integration by parts and Girsanov's 

theorem, taking Q*'^-conditional expectations, we have Sf {OVt^ + {I - e)Vi^ - Vf) < 0, 
which gives the result. □ 



4 The second optimization problem 

In this section, we assume that Us = U (c^) and Ut = U (ip) where U and U are given utility 
functions, c is a non-negative G-adapted process and ip a ^^-irieasurable non-negative 
random variable. We fix a probability P equivalent to P with a Radon-Nikodym density Z 
with respect to P given by: 



n 

dZt = Zt^{etdM^ + 5^(e-^* - l)dNi), Zo = 1 
1=1 



(25) 
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4.1 Formulation of the problem 

We study the following the optimization problem of the robust maximization initial prob- 
lem (1): 



sup E*- 

{c,V)G^(x) 



+ E'^ 



(5,5f In. 



ds + Si In Z, 



:= sup Vq 

where A{x) is a set of admissible positive control parameters, Vb is the value at initial 
time of the value process V, part of the solution {V, M^''^, v) of the BSDE (6), in the case 
Us = U{cs) and Ut = tJW- Here, Q* is the optimal measure for V{U{c),U{'il:)), and 
depends on (c, The preferences are modeled by the utility functions U and U which 
shall satisfy the following conditions: 

Assumption A 3. The utility functions U and U defined on IR+ satisfy the usual conditions: 

i) Strictly increasing and concave. 

ii) Continuous differentiable on the set {U > — oo} and {U > — oo}, respectively. 
Hi) U'{oo) := lim^^^oo U'{x) = and U'{oo) := lim^^^oo U'{x) = 0, 

iv) U'{0) := lima;->o U'{x) = +oo and U'{0) := lima;->o U'{x) = +oo, 

v) Asymptotic elasticity AE{U) := lira sup < 1. 



The set A{x) of admissible control parameters is given by : 



Definition 5. 

A{x) is the set of control parameters (c, V') G 'H^([0, T]) x L^(il, Gt) such that 



E^[ / ctdt + ^p] <x, (26) 
Jo 

and {U{c),tJ{ip)) G Z^f" x L^^p and {cU'{c),iptJ'{Tp)) e x L'^^'p for any pair {c,tp) G 

T-L'^{[0,T])xL'^{Q,QT),asivellastheprocessexp{'-f \U {ct)\dt) (respectively exp (7 |ct||C/'(ct)|(it)) 
belongs to the class [D] (see Dellacherie and Meyer, pp.89. Chapter VI [12] for definition). 

Remark 6. Under our assumptions, we have the set A{x) is convex and closed in the topology of 
convergence in measure (see, Cuoco [11] Lemma B3., pp.70). 

In order to clarify and motivate the constraint (26) satisfied by the control parameters 
(c, tf^), we present a generic example in a financial market where the process c can be inter- 
preted as a consumption and ^ as a terminal wealth : 

Example 2 (Consumption-investment problem). We assume the same model as in Example 1, and 
we consider a financial market consisting ofd + p+1 assets. The savings account is assumed to be 
contant equal to 1, the prices of the d + p risky assets are G-semi-martingales given by 



dSl = SI 



li\dt + 



d 

E 



dNf 



k=l 



t 



I,..., d + p 



(27) 



Robust utility maximization problem 



17 



where N = {N^,. . . ,N'^) is the martingale Nl = l^.<t - /g^^' \\ds and W = {W^, . . . , WP) 
is a Brownian motion; a is the {d + p) x d volatility matrix {a^''',i = I, . . . ,d + p;k = 1, . . . ,p) 
and if is the (d + p) x d matrix {f^'^ ,i = 1, . . . ,d + p;j = 1, . . . ,d). We note also note S the 
(d + p) X (d + p) matrix S = [a, tfX]. 

Given an initial wealth x, the wealth process {X^'^''")Q<t<T associated to the triple {x, c, vr) where 
vr zs the portfolio strategy and c the consumption plan, follows the dynamics given by: 

dX^^^'"" = TrtdSt-Ctdt, X^'^'^ = x, (28) 

The set of consumption-investment strategies (c, vr) satisfying the following no-bankruptcy condi- 
tion is called the admissible strategies set and denoted by A{x) : 

P - a.s., Xf'"'^ > 0, Vt e [0,T]. (29) 

For this model, one assume that : 

• The appreciation rates {fi^ ,i = 1, . . . ,d + p) are bounded predictable processes. 

• The processes {f^'^ , 1 < i < d + p,l < j < d) are bounded and predictable and satisfy 
ifl'-' > — 1 a.s. 

• The processes (cr*''^, 1 < i < d + p, I < j < d) are bounded and predictable 

• The matrix S is bounded and invertible. This condition ensures that the market is arbitrage 
free and complete. 

Then, the pair consumption-terminal wealth satisfies the budget constraint 

r ctdt + X'/^^] <x 
Jo 

where P zs the unique equivalent martingale measure for the state price S with density Z given by 
(25), where [6*, e"^ - 1] = S" V • 

4.2 Properties of the value process 

In this section, we derive regularity properties for the value process V^''^^'^ with respect to 

Proposition 4. Define the map G : A{x) — > Dq^^ as G{c,ip) = V, where {V,M^''^,v) is the 
solution of the BSDE (6) associated with {U{c), U{tp)). Then 

(i) G is concave, i.e., for all 9 G (0, 1) and (c^,^^), (c^, V'^) € A{x): 

G {9c^ + (1 - 0)c2, eij^ + (1 - 0)7/;2) > eG{c\ij^) + (1 - 9)G{c^,ij^). 

(ii) Let Go{c, ip) be the value ofG{c, ip) at time 0, i.e., Go{c, ip) = Vq. If{c"-,ip'^) G A{x) converges 
decreasingly to (c, ip) G A{x), then Go{c'^,^l)''^) converges decreasingly to Gq{c, ip). Moreover Go 
is upper continuous with respect to the control parameters. 
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Proof: Let M 



be the solution of the BSDE (6) associated with (C/(c'=), U{i;'')) 



for k = 1,2. For any 9 G (0, 1), let (F^, M^''^, u^) be the solution of (6) associated with 



iUi9c^ + (1 



, C/(6'V'^ + (1 - and (F'', M^'^ v^) be the solution of (6) associated 



with 9U{c^) + (1 - 9)U{c^), 9U{ij^) + (1 - 9)U{^p^)) and set F := 9V^ + {1 - 9)V^. Then, by 
using both the concavity properties of (U, U) and Theorem 2, we get >V^. Moreover, 
as consequence of Theorem 3, we obtain > V, which gives the assertion (i). 
Let us now consider (c",V'") a decreasing sequence of control parameters in A{x) such 



ct, Vt a.s and -0" — > ip a.s; then, by using inequality (19), and the fact that the 



that 

functions U and U are non-decreasing, we get 

rT 



•^0 I — 



{uid:)-uics))ds + {u{r)-um 



(30) 



where Q* is the optimal measure associated with {U{c), U (ip)). Thus, by using the mono- 
tone convergence theorem V'^"'^" converges decreasingly to V^'^. Let (c",V'") € A{x) 
he a sequence of control parameters such that c" — > c a.s and V" — > a.s where 
(c, '0) G A{x) and denote c" = sup^>„ c"^ = sup^>„ Then, — > c a.s decreas- 
ingly and ip^ - 
and therefore: 



a.s decreasingly. It follows that converges to V^''^ decreasingly 



< lim^n ''^ 

n 



^0 



limsupvQ 

n 

Hence, Go is upper semicontinuous with respect to the control parameters. □ 

Definition 6. The pairs {c^,%l:^), {c^,il)'^) G A{x) are comparable if either (c^, i/;^) > (c^, ■0^) or 
(c^, lA^) < (c^, V'^) wiih the order introduced in Definition 4. 

Proposition 5. Assume that Assumption A. 3 holds and let (c^, -0^), (c^, V'^) be two comparable 
plans in A{x). Then the function defined on (0, 1) and valued in D^'^ 



is right continuous at 0. 



Proof: Assume first that (c\ V'^) < (c^ V^)- Let, for e g]0, 1[, = G(c^ + e{c^ - c^), + 
e(V'^ — V'^)) and F = G(c^,V'^)- From Proposition 3 and the obvious inequalities U{c^ + 
e(c2 - c^)) > C/(ci) and tj{'^^ + e(02 _ ^i)] > [/(^i), we obtain 



sup \Vt 

0<t<T 



< CK 

-T 



r.2 



+ / \Uicl + eicl-cl))-Uicl)\'ds 







Using now the concavity properties of U and U, we obtain 

< U{cl + e{c',-cl))-U{cl)<eU'{cl){c',-cl) 

Thus, we have 



E^ 



sup 

'0<t<T 



Vt - Vf 



< GE^ 
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Assume now that (c^, V'^) > (c^, ^^)- Then, using the fact that G is concave with respect to 
the control parameters, one has > (l — e)y^ + ey^ where are associated with {c^,ip^), 
hence > _ Moreover, since + e(c^ — c^) < and -0^ + e('(/'^ — ^"^) < 



we have by Theorem 2 that > 
now Proposition 3, we get 



E^*''l sup 

'0<t<T 



. Therefore 



< \Vt^-Vt^\,t G [0,T]. Using 



supo<t<r 



Finally, we conclude there exists a constant C > such that: E''^* 

where Q* = Q*-^ if {c\ip^) > (c^,^'^) and Q* = Q*'^ if (c\V^) < (c^V^), then by Kol- 
mogorov's criteria, we deduce that ^' is right-continuous at 0. □ 



We give now a differentiability regularity result for our BSDE's (6) . This result is not well 
known even for a standard quadratic BSDE's. The proof will be given in the Appendix. 

Theorem 4. Let (c^, and {0^,11?) he two comparable plans in A{x). Let (V^, ■u'^) be the 
solution of (6) associated with {U{c^ + e{c^ - c^)), U{ip^ + e(V'^ - and M^'^, v^) the 
solution of (6) associated with (C/(c^), JJ{iIj^)). Then, is right-differentiable with respect to e at 

0. Moreover, if we denote by d^V := lim 



-, then there exists d^M^'^, d^v G 



CW Q^'*) such that the triple {d^V, d^M^'", d^v) is the solution of the following BSDE: 



dd.Vt = {5td,Vt - U'{c\){ci - c\)) dt + dd,M^^' + devldNl C 

i=l 

[d,VT = U'{i;'){i;'-i;'), 
where A* := A*e~^'^'' and TV* := TV* - /ole""' ' - l)Xldt is a Q^^* -martingale. 
Moreover, we obtain for all t < T: 



a.s. 



(31) 



d.Vt = E^ 



(32) 



4.3 The Dynamic maximum principle 

We recall that we are interested in the following optimization problem: we associate with 



a pair (c, -0) G -4.(x) the quantity X^ 







E* 



/q Csds + and we study 



u{x) = sup 



(33) 



Here V; 







Vb, where (F, M^^ u) is the solution of the BSDE (6) associated with {U{c) , U{ij)). 



Note that if we are in the setting of Example 2, our problem correspond to a maximization 
of recursive utility function over consumption-investment strategy where Xq is the initial 
wealth associated with the consumption c and terminal wealth ip. 
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Proposition 6. There exists an optimal pair {c^,ip^) which solves (33). 

Proof: The uniqueness is a consequence of the strictly concavity property of Vq. We shall 
prove the existence by using Komlos theorem. 

First step: Let us first prove that sup(c^)g^(^) Vq'^ < +oo. Because F G Qj, we have: 

r- r'^ 1 

sup Vq'"^' < sup E'^ U{i;) + / U{cs)ds := u{x) 

Using the elasticity assumption on U and U, we can find 7 e (0,1) and xq G K such that, 
for any 9 > 1, one has: 

u{ex) < e'<u{x) Vx > xo, 
u{ex) < e^u{x) vx > XQ, 

hence, for any x > xq: 

u{ex) = u{e-^) + / u{e^)ds < e^x)- 

where (c^^W^^) is the optimal control of the optimization problem under the set of mea- 
sures Q = {P}. Hence, AE{u) < 1, which permits to conclude that, for any x > , 
u{x) < +00 (see [?] and [38] chap. 3, Lemma 3). 

Second step: Let (c", 4^"^) € ^(x) be a maximizing sequence such that: 

lim Vf'^" = sup V^'''' < +00, 

where the RHS is finite thanks to step 1. Using Komlos criterion, we can find a convex 
combination (c", '(/'"■) G cont;^(c"', ^"), (c""'"^, V'"^"^), ■ ■ ■ ^ which converges P-a.s. We de- 
note by (c*, V'*) this limit, which belongs to ^(x) since it is a closed convex set. Moreover, 
there exists Nn > n and a positive sequence (^™)mgN satisfying Yl,rn=n ^"^ = 1 such that 
{cn,'4'n) = (X]m=n ^'"c'^i Z]m=n Therefore, by using the concavity and the increas- 

ing properties of the functional Vq which respect to the control plan we get: 

,i,n \^Nn nm m y^iVn ffm jm . ^ rn ,/,m _n ,/,n 

'''' = Vq'^=" '^m=n '/^ > 6* > 

r?i=n 

Moreover, using the upper semi-continuous property of the functional Vq which respect to 
the control plan we get: 

sup Vq =limsup\/Q <limsup\/Q = Vq ^ . 

{c,ip)GA{x) n n 

□ 

In order to characterize the optimal solution, we recall the classical convex analysis result. 
Proposition 7. There exists a constant u* > such that : 

u{x) = sup I Vi^'"^^ + u*(x- xi^^'^A I (34) 
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and if the maximum is attained in (33) by {c*,ip*), then it is attained in (34) by {c*,ip*) with 



X, 







X. Conversely, if there exists > and {c^,ip^) such that the maximum is attained 



in 



sup 



^ + z/° ( X - X 



with Xq = X, then the maximum is attained in (34) by (c , )• 

Let v > Ohe fixed and L be the map given by L(c, tp) = Vq'^''^^ — uX^^'^K We now 
study the following optimization problem: 



sup L(c, ip) . 



(35) 



Proposition 8. The optimal plan {c^, il^^) which solves (35) satisfies the following (implicit) equa- 
tions: 



7^ ij 



U'{d 



, dt®d¥a.s, U'{^p'') = ^^,dFa.s 



(36) 



where is the Radon-Nikodym density of the probability measure associated with the optimal 
plan (c°, V'^). 

Proof: Consider the optimal plan {c^, which solves (35) and another plan (c, ijj). For 
e G (0, 1), one has L(c° + e(c - c°), ^° + e(^ - V°)) < ^(c°, then 



(37) 



From the definition, we obtain that 



d^xf'^'^ :=lim-(Xo 



1 - ,.(cO+e(c-cO),V''+e(c-cO)) _ ^(cO,,/<0)^ 



(c, - c'Dds + (^ - V°) 



Taking the limit when e goes to in (37), we obtain: 



< 



(38) 



where deV^^°'^°'' exists and is given explicitly by Theorem 4. From the explicit expression 
of we get: 



T 



S'TZ)^U'm{i,-^^)+ / 5,^Z,"C/'(c^)(c.-c^)ds 

Jo 

Jo 



It follows from equality (38) that 



S'^z!^U'{^lJ')-uZ^]ii,-r)+ 1^ {StZ^U'{c^,)-uZ^]{cs-c^Ms 



T 



< 
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The end of the proof is the same as in El Karoui et al. [18] (proof of Theorem 4.2, p. 677). 



In particular, for any i/j, ¥F 



S^Z^U'iip^) - = a.s 



< 0, hence 



We find the optimal c with similar arguments. 



□ 



Theorem 5. Let I and I be the inverse of the functions U' and U'. The optimal plan {c'^,ip'^) 
which solve the problem (34) is given by: 



,0 vP 



Z 







Z 



c^ = /(_-^) dt^dFa.s, ip'' = I \ ] a.s 



q5 yO 

iJrp ZJrp 



where > Q satisfies: 



T I ,,0 



v'^ Zf 

sfzf 



dt + I 



X. 



Proof: Define the map: / : (0, +oo) — > (0, +oo) as 



T 



s{ z? 



dt + I 



c5 yO 



Then, using assumption A.S, / is monotone and satisfies lim,^_i.o f{^) = +00 and Imiy^j^ao /(z^) 
0. For any initial wealth x € (0, +00), there exists a unique such that /( 



V ] = X. 



Let (c, V) e A{x) and {V^^'^'> , M^'^ v) (resp. v°)) the solution of the BSDE 

(6) associated with (U (c^), U (^°)) (resp. {U{c), U {ip))) then from the inequality (19) (see 
the comparison theorem), we get: 



S't (C/(V) - U{^')) + / 5f {U{cs) - C/(cO)) ds 



S'TU'{r){i^ - O + / 5f f/'(c^)(c, - c",)ds 



&TTl(J)\ 



It follows that: 



(y¥ f-T yF \ 



T 



Since (c, V') G A{x), then E^ ip + Csds < x. Using that E"^ V° + Jo c°fis 
conclude: 

^0 — ^0 



X, we 



□ 
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5 Logarithm Case 

In this section, we assume that the process 6 is deterministic and that U (x) = ln{x) and 
U{x) = 0, hence /(x) = — for all x S (0, +oo). We introduce, as in Theorem 5, the optimal 



process Cj = H—-^) = — - — • Recall that the Radon-Nikodym density Z, and the 

\S^ Z^ J ^ Zt 

Radon-Nikodym density of the optimal probability measure Z* (given in (9)) satisfy 

n 

dZt = Zt^iOtdM^ + - ^)dNl), Z^ = 1 (39) 

i=l 
n 

dZl = ZU-dM^'' + Y,{e-''''-'^)dNi),Z*o = l. (40) 

1=1 

For any deterministic function a such that a(T) = 0,V admits a decomposition as 

Vt = a{t)\n{c;) + pt 

where /3 is a process such that /3t = 0, for a specific choice of a. Our goal is to characterize 
the process /3. As in [6], we introduce Jt = i:^^^f^t in order to obtain a simple BSDE. Note 
that, even if Z* is implicit (the coefficients depend on the solution c*), the BSDE for J is 
explicitly determined in terms of the given parameters A* and of the given probability P. 

Proposition 9. The value function V has the form Vt = a{t) In(cJ') + (1 + a{t))Jt where a{t) = 
— e-A and {J,M-^'^,j) is the unique solution of the following Backward Stochastic 

Differential Equation: 

dJt = ( (1 + + Kt))Jt - k{tW)) dt + dM^'^ + '^d{M^^^)t + \k{t){l + k{t))0U{M' 

d d 

+ Y.AdNl + {g{jl)K + {m{e-'^ " 1) + e""^'^^^ - l) Aj)di 

i=l i=l 

Jt = 

where k{t) = — j^^^. Here, the processes M^'^ and dNl = dHl — \\dt are F-martingales where 

dF\g^ = ZtdF\g„ Xi = e'^W^Uj and 

dZt = -Z,-(^{t)etdMf-Y,y^'^^l-l)dNl^ (41) 

Note that, in a complete market, one obtains a forward backward system for the pair J- 
optimal wealth. 

Proof: Using the fact that V satisfies the BSDE (6) and the assumed form of V in terms of 
(a,/3), one obtains dVt = {5{t)Vt -\n{c*t)) dt - d{\D. Z"^) = a{t)d{\ndl) + {hicl)a' {t)dt + dPf 
Therefore 

dpt = S{t){Vt + a{t))dt - (1 + a {t))\i\{ct)dt + a{t)dhiZt + (a(t) + l)dlnZ; 

= ( {5{t)a{t) - 1 - a'{t)) lnc*(i) + 5(t)/3f + a{t)S{t))dt + a{t)dln % + (a{t) + l)dlnZ; 
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We choose a so that 6{t)a{t) = l + a'{t). It follows that 

d/3t = d{t){f3t + a{t))dt + a{t)dlnZt + (a(t) + l)dlnZ; 
After some obvious computations taking into account the form of Z and Z*, one obtains 

d 

d(]t = 6{t){Pt + a{t))dt + ^(^{a{t) + l){e-yi -l)-a{t){e-''^ -l)^Xidt 

i=l 

+a{t)9tdM^ + (a(t) + l)dMf'" - ]^[a{t)e^^d{Nnt - {a{t) + l)d(M^''^)t) 

+ ^((a(t) + l)y,'-a(t)4)dF^ 
j=i 

We now define J = jTp^^/?* and set k{t) = — j^^^, then we find the following dynamics: 

+dMj^''= - k{t)etdM^ + - {k{t)eld{M^)t + (i(M^'^)t) + ^ (yj + k{t)zl)dNi 

i=l 

We introduce the martingale M'^''^ as dM/''' = dM^^''' - k{t)dtdM^. It is easy to check that 

d{M'^^^)t = d{M^'^)t - k^{t)9^d{M^)t - 2k{t)etd{M-^'^,M^)t 
and we denote jl = y't + k{t)zl. Using the fact that, due to the form of g, for any x, k, z, X, 
xdN + X{g{x-kz) + kg{z))dt = x{dN - {e''^ -l)Xdt) + (^g{x)e''^ + (e"^ - l)k + e''' - l) Xdt 
one obtains 

d 

dJt = ((1 + 5{t)){l + k{t))Jt - 5{t)k{t)yt + (5(jl)e'=(*)"' + k{t){e-'' - 1) + e^^^^' - l) Ajdt 

i=l 



+dM^^^' + ^d(M-^'=)f + k{t)etd{M^''', M'')t + \k{t){k{t) + l)0t2d(M'=)f 



i=l 

k{t)etdM^ - Y.t=M^^^'' - '^)dNl\ . The 



We define P as = Zd¥, where dZt = -Zt 

processes M-^'^ and defined as dM^''' = dM^'"" + k{t)etd{M-^''', W)t and diY^^ = dNf - 
_ I'^y-^dt = dHl — X\dt are P martingales. The result follows. □ 

6 Appendix : Proof of Theorem 4 

Let M"'^, v^) be the solution of (6) associated with {U{c^+e{c^-c^)), U {ip^ +e{ilP -^^))) 
and (y^, M^'"^, t;^) be the solution of (6) associated with [U (c^), ^/(V'^)) and denote 

:= , A.M-^ := , A,?;* := , 

^ ^ _ ^ _ (42) 
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then, (A^y, A^Af^, A^v) satisfies the following equation: 



A,Vt - [\dsA,Vs - A,Us)ds = ^{{M^'^h - {M^'')t) 
Jo 

+ -E [\9ivr)-9{vln)Kds + A,M,' + Yl f ^^<dNl 



(43) 



with final condition A^Vt = AJJt- 

We first start by giving the following a priori estimates for the sequence {A^V, A^M'^, A^v): 

Lemma 1. Assume the same conditions as in Theorem 4. Then, there exists a constant C > such 
that: Vi = I,-- - ,dyp e N*,Ve > 0, 



sup 

0<t<T 



\a,v,U{a,m^)^ + Y.£ ■ 



A,vi 



-Xlds 



<C, 



(44) 



where A^M^ is the Q*'^ martingale part of the Q*'^ semimartingale A^M*^, and A* := A*e is 
the intensity process of the process W under the probability measure Q*'^. 

Proof: Let (c^, V'^) and (c^, V'^) be two comparable plans. We introduce the processes 



Kl := 



exp 



{5sV! - U{c\ + e{cl - c\)))ds - + e{^^ - i^^)) 



Kl := E^[exp^ {5sV}-U{cl))ds-U{i;'] 

Obviously, for all t G [0, T], one has: 

= -ln{K^)+ [\6,V,^-U{cl + e{cl-cl))ds 
Jo 



HKl)+ f{6sV}-U{cl))ds, 
Jo 



hence. 



In 







W] 


Jo . 



- {V: - V}) - - {U{cl + e{cl - c\)) - U{c\)) 



ds . 



(45) 



K 



For t e [0, T], we define ■= and Kl = [Kn . The processes K^ and {K^)'^ are 

Kf V / 

positive semi-martingales which belong to LP(P) since: 

{Kl)-P = exp + fp{U{cl + e(c2 - a])) - SsV:)ds . 
Jo 
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In the other hand, by using the dynamics of K"^ and under the probability measure P: 



dKt = Kl_ 
dKl = Kl 



-dM:'' + ^{e-<" -l)dNi 

d 

-dMl^' + ^{e-^'-' -l)dNl 



i=l 



and applying integration by parts formula, we get the dynamics of given by: 

dK^ = Kf_[- d{Mt''' - mI'^ - {M''^ - M^'^, M^'^)t) 

+ (e-(^*'"-^''') - 1) [dHi - e-^^'Ajdi] ] 



Clearly, K'^ is Q*'^-local martingale. Then, the processes K'^ and {K*^) ^ are positive ' 
submartingales. We now split the study into two cases. 

First case: (c\ V'^) < (c^, ^^). Using the inequality (19), for all t G [0, T]: 



(46) 



\^eVt\ < E' 

1 V 



T qS 



- V^^) + / ^U'{ct){cl - cl)ds 



t sf 



sup V 
0<t<T \J^t 



< exp 



p(||<5||oo + l) sup |AeVt|+ r pU'{c]){cl-c\)ds 

0<t<T Jo 



Setting K = p(||5||oo + 1)/ we obtain: 



sup (wiY ^ f'^ sup E^"" [C7'(^i)(^2 _ ^1) ^ U'icDici - cl)ds 

^<t<T \J^tJ ^ 0<t<T Jt 

+ j\u'{c]){cl-c\)ds 



Using Jensen inequality, we have: 

sup { ^r] ^ sup E"^ 

o<t<T \ Kl J 0<t<T 



exp U'ii^'){r -r)+ U'{cl){ci - ci)ds 



xeMf pU'{c]){cl-c])ds). 
Jo 



(47) 



Thanks to the assumption (c*, -(/'*) G A{x), we conclude that supg^^^-^ \^) ^ '^^^^)- 
Second case: (c^, V'^) > (c^, V'^)- Then, using concavity property, we obtain for all t G [0, T]: 

^.Vt\ < \V,' - V,% \A,U{cj)\ < U'{c^,){cj - c?) 

Now, using the same arguments as in the first step, we get that: 



sup 



0<t<T \Kl 



< sup E^ ' 



0<t<T 



exp ( U'{^^){^' -^'')+ I U'{ci){ci - ci)ds 



T 



X exp( / pU'{ci){ci - ci)ds) 
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We use the same arguments to prove supo<«'j- \K^\ G LP(P). From the representation 
theorem, there exists two continuous martingales M'^ '^, M'^''^ and d predictable processes 
, A;^ such that the martingale parts of and ~ are given by : 

i=l ^0 i=i ^0 



kldNl 



These processes and ~ being positive Q*'^-submartingales, using (47) there exists two 
constants Ck and Ck such that: 



{klYKds 



2" 



1 



T 



(48) 



<Ck 



From the uniqueness of the representation theorem and equation (45), we get, for 1 < i < 



1 + 



and /^tvl = In 



d, —^evl = In 
-^T^ + \k^ \K^_ . Moreover we find: 



1 + kiK^. .Therefore we find exp(|A,7;|j) - 1 < 



< 







^-^Kds+ I \kl'^\KIXlds 







1 



t<T Jo 

Using Cauchy Schwartz inequality, we find that: 



sup^ / |fc^'*|A*ds + sup-fCM \kl'^\Xlds 



t<T Jo 



E^ 



< 



+ 



E«3*'Msup. 



1 



^t<T {Kir Jo 



T 



E^' sup(Kr: 

V t<T 



Xlds 



Xlds 



E*^"^ /^IM'f Aids 







E^ 



Let prove now X\dt is square integrable under the probability Q*'^. We write first the 
expression under P using Bayes formula: 



E^ 



T 



Xldt 



Z. 



T 



e-"^''Xlds 







<E^ 



z: 



T 



T 



T 



Xlds I e-^^'^'^Js 



then we use Cauchy Schwartz inequality to find the following estimates: 

rT 



E^ 







r Kdt 


<Cl( 


Jo 





E 



e-^''''''Kds 







where we make use several times of Assumption A2-iii). Moreover, we can see that 



E 



= E 









- 1 + 1) Xlds 



< 16E 



T 



1) X\ds+ / X\ds 
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Therefore, since the martingale -M^'"" + / J^f^^ [e-'"^'' - 1 j dN} belongs to LP(P), and by 

assumption A2-iii) again, we conclude that E'^ Jq {e~'"'' — 1) A* ds < +oo for any p > 1. 

Moreover since Z^*' e Lp(P), we get that E^*' A* ds < oo. Then using again Cauchy 
inequality: 



< C E 



1 

sup — J 

0<t<T {KD 



2 /_,0*"1 



\K''\ Kds 



+C K"^ - 



sup (i?,^)^ 

0<t<T 



2 



|A;f I Aids 







From (47) and (48), we deduce that there exists a constant C2 > such that: 

rT 



E*^ ' 



l)\lds 



and then using the expansion of the functional x — > we get: 



E^ 







< C2P\. 



(49) 



In order to conclude the proof of the lemma, it remains to establish that there exists a 
constant Ci satisfying: 

E^*''[(A,M^)t] < Ci. 



First case: (c^,^^) > [c^ ,i)^),\henU{c^+e{c^-c^)) > [/(c^) and [/(^i+e(V'^-V'^)) > Uitp^). 
From Proposition 3, it follows that: 



E^ 



< E*^ 



0<t<T 



i=l 



sup \Vf - Vl\' + (il/^'^ - M''^)t + I {vr - vVfy^tdt 



Since 
and 
we get: 



< U{cl + 6(c? - c^)) - U{c\) < eU'{c\){cl - c\) 



sup |A,T4|2 + (A,iVr)T+ / {A.vlfxids 

0<t<T 



< E^ ' 



T 



(i7'(^i))'(V'2-V^)' 



+ / {U\cl))\c',-clfds 
Jo 
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The process Z'Q*'' belongs to LP(P);moreoverC/'(V'^)(V'^-V^) e L'^^'p and U'{c^){c^ -c^) £ 
D'^^ since (c^, ^^), (c^, G A{x). It follows that there exists a constant C > such that: 



sup jAeVtj 
0<i<T 



< C 



Second case: (c^-iA^) < (c\V^)- We first prove that for all t G [0,r], > 1. Let us 
recall that: 



= exp {^-A,Vt + i5sA,Vs - A,Us)ds 
Define the process X as 

Xt = -A,Vt + [ {6sA,Vs - A,Us)ds, 0<t<T. 
Jo 

From integration by part formula, we get: 



S^Xt 



-A,Vo 



StdA.Vs 



StAM.ds 



Since the process 6 is positive and bounded, there exists a constant L > such that < 
L <1. It follows that: 

SfXt > (-1 + L)A,Vo - LA.Vt - [ SiA.Usds 

Jo 

Note that, for all t G [0,T], A^Ut < since (c^,^^) < (c^,^^) and using comparison 

theorem A^Vt < 0. Therefore, for all t G [0, T],Xt> 0. Finally, > 1. 

In the second step of the proof, we give the dynamics of the process K"^ using Ito's calculus: 

dKt = Kl_ -dA,M^ + ^(e — - l)dNl + dAt 



i=l 



where A is an increasing process. Since K'^ is a positive Q*'^-submartingale, we obtain 
from (47) and Kl>l: 



(A,M" 



I J. 



then we conclude: 



E^ 



{KlYd{A,K'nt 



< E^'- 



<Ck 



(A,M^)t 



<Ck. 



Finally, by using concavity property we have shown that: lA^Vf | < 1^4^ — for all t G 
[0, T], then: 



E<0*'' 


sup |A,Vt|2 




sup 


< 2EQ*'' 


1 T 1 2 1 T r2 1 2 

sup l/j + sup l/j 




te[o,T] 




te[o,T] 




tG[0,T] te[o,T] 



Therefore, since the process usmg 
Cauchy Schwarz inequality that there exists a constant C such that: 



E^ 



sup \AA^t\ 

iG[0,T] 



< C. 



□ 
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Proof of Theorem 4: Let recall first the equality : 
then the equation (43) may be written as: 

= J_(M^.c _ ^i,c^^ + V r [-(e-^"' - e-'^'-') + e-'^''' \\ds 

+ (A,Mf + (A,M^ Mi'^)i) + V r A^vlidNi - {e-^'" - l)\\ds). 

ZTT Jo 



/(^(e — l)X\ds are Q^ *— martingales. It follows that the process: 



By Girsanov theorem, the processes := A.M" + {A.M" , M'^'") and A^^ := A^* 

rt 

'0 t>0 

is a Q^'*-submartingale and we have the following decomposition: 



{A,Vt - f{5sA,Vs - AJJ,)ds) 
Jo t 



A,Vt - [\dsA,Vs - A,Us)ds = UA,K-nt + f [-(e"^"' - e""' 
Jo ^ .^1^0 e 

+ e^'^'^A^vi] Xids + AeMf + V / A.v^.dNl 

i=i Jo 



(50) 



Moreover, using the uniform estimate (44), we get: 



linjE«-'(l(A,AnT)<Climi = 0, (51) 



and using the expansion of the functional x — )■ e^, we get from (49): 



thus, passing to the limit as e ^ 0, we conclude that: 



limE'3*'' ( /"^[l(e^"^'' - e-'''"') + e"^'" A.-y^JA^ds ) = 0, 1 < i < d. (52) 
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Moreover, the estimate (44) ensures that the sequence (A^V, A^M'^, A^v)e>o is bounded in 
?^^(P)x7Wq(P) x£2(A, P). As a consequence, we can extract a subsequence {A^^^V, A^^.M'^, A^j^v)^^^^ 
which converges weakly in W^(P) x 7Wq(P) x >C^(A,P) and by Banach-Mazur Lemma, 
one may construct a sequence (V^^, M'^''^, ?}^)e>o of convex combinations of elements in 
(A,,y, Ae,M^, A,^v)keN of the form 



j=\ j=\ j=i 



such that M^'^ v^)^>o converges strongly in ?^2(p) x A4^(P) x/:2(A, P) to [d^V, 5eM^ d^v). 
Moreover, the triple (F'^, M*^'*^, 9*) satisfies the BSDE's (50) associated with {V, U ) where 

i=i j=i 

Therefore, passing to the limit in this equation, thanks to (51), (52) and the dominated 
convergence theorem, we get that [dfV^ d^M'^, d^v) solves the BSDE's 



dd.Vt = {hd.Vt - U'{c\){cl - c\))dt + dd,M^ + devidNl O.Vt = U'O^'Xi^^ - i^^)- 
Therefore (Sfd.Vt + Jo S^,U'{cl){c1 - cl)ds)^^^ is a Q*'i martingale which can be written 



as: 



Std,Vt+ I 5f[/'(cJ)(c2 -cj)ds = E^^' S'Td,VT+ I S',U'{c',){ci - c',)ds 



Gt 



Hence we get: 



deVt 



T 



o5 



□ 
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